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We study a model of "elastic" lattice polymer in which a fixed number of monomers m is hosted by 
a self-avoiding walk with fluctuating length I. We show that the stored length density pm = 1 — (0 /m 
scales asymptotically for large m as pm ~ /5oo(l — 6/m + . . .), where 6 is the polymer entropic 
exponent, so that 6 can be determined from the analysis of pm- We perform simulations for elastic 
lattice polymer loops with various sizes and knots, in which we measure pm- The resulting estimates 
support the hypothesis that the exponent 9 is determined only by the number of prime knots and 
not by their type. However, if knots are present, we observe strong corrections to scaling, which 
help to understand how an entropic competition between knots is affected by the finite length of 
the chain. 

PACS numbers; 02.10.Kn, 36.20.Ey, 36.20.-r, 87.15.A- 



I. INTRODUCTION 

According to renormalization group theory, the scal- 
ing properties of critical systems are insensitive to micro- 
scopic details and are governed by a small set of universal 
exponents |1] . Also polymers can be considered as critical 
systems in the limit where their length I (the number of 
chained monomers) diverges PHI]. For instance, the ra- 
dius of gyration of an isolated polymer in a swollen phase 
scales as Rg ^ l" , where z/ « 0.587597(7) [5] in d = 3 di- 
mensions is a universal critical exponent. One of the sim- 
plest models in the universality class of swollen polymers 
is that of self-avoiding walks (SAWs) on a lattice. Hence, 
these have been used extensively to extract information 
on critical exponents and scaling functions ^2 , 4- 23J . The 
total number of SAWs, i.e. their partition function, has 
the following large-Z expansion 

Zi^^iH'{l + Al-^ + ...). (1) 

Here non-universal (model-dependent) quantities are the 
connectivity constant /i and the amplitude of the cor- 
rections to scaling A. The entropic exponent 9 de- 
pends only on boundary conditions: in d = 3 we have 
6 = -f — 1 = 0.1573(2) [24J for an open chain whereas 
9 = a-2 ^ -dv = -1.762791(21) for self-avoiding 
polygons (SAPs), that is, linear chains with the two ends 
on adjacent lattice sites. Renormalization group analysis 
suggests that the exponent A, characterizing the leading 
corrections to the scaling behavior, is also universal ^j |3] 
Models with full self-avoidance, such as SAPs, have 
been used to study the statistical properties of knotted 
chains ^13] [T4j [T9l [251 - 135] . Knots in polymers have at- 
tracted a lot of attention during the past years, also be- 
cause of their occurrence in biopolymers as DNAs, RNAs 
and proteins [36-41 . As usual, SAWs represent a mini- 
mal effective model to grasp the essential, coarse-grained 
features of polymer chains. Simulations of knotted SAPs 



in ensembles with fixed topology are performed with a 
grand-canonical algorithm (BFACF [6 , from the name of 
the authors) tuned to span a range of chain lengths (al- 
gorithms with fixed N are not ergodic in this case). For 
this algorithm, the tuning of step fugacities to urn l//i is 
necessary to achieve samplings of long chains. It would 
be desirable to have a simpler and more stable method 
to sample the same chain lengths. 

In this paper we study a class of polymers referred to 
as the elastic lattice polymers (ELPs), which are SAWs 
accumulating some stored length along their contour. 
This leads in fact to a partial lifting of the self-avoidance 
condition between consecutive monomers of an ELP, as 
sketched in Fig. [l] We will consider equilibrium prop- 
erties of polymers with a fixed number of monomers m, 
and in which as a consequence the length ^ < m of the 
self-avoiding backbone described by the monomers fiuc- 
tuates. This explains the name "elastic", and implies a 
resemblance with the class of grand-canonical SAW mod- 
els. 
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FIG. 1. Example of an elastic lattice polymer on a square 
lattice. This polymer is composed by m = 11 monomers 
describing a SAW backbone of length I = 8 (dashed area). 
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There are several reasons for studying this model. On 
the theoretical side, it can be considered as an enhanced 
SAW: besides sharing critical exponents with SAWs, its 
fluctuating length enables new avenues to estimate criti- 
cal exponents. ELPs have been used in studies of polymer 
dynamics as phase separation in polymer melts [22_, or 
in translocation through nanopores [21 . but their equi- 
librium properties have so far received little attention. 

The key quantity we focus on is the equilibrium aver- 
aged stored length density defined as 



(0 



(2) 



where (Z) depends on m. As will be shown. has 
a simple asymptotic behavior for large m from which 
one can extract universal exponents: the leading cor- 
rection to the asymptotic value for p,„ scales as 6/m, 
where 9 is the entropic exponent defined by Eq. ([T]). 
We illustrate the result of this approach for the case of 
ELPs with fixed knots. If knots are present, the stored 
length approaches its asymptotic value with strong, knot- 
dependent corrections to scaling. The expectation of a 
homogeneous stored length within an equilibrated chain, 
combined with the knowledge on how its density varies 
with the chain length, leads us to a new view on the issue 
of entropic competition of knotted regions [23 • On the 
numerical side, we find that ELPs, compared to grand- 
canonical algorithms, have the nice feature of stabilizing 
the sampling quite narrowly around an easily tunable 
length (l). 

This paper is organized as follows. In Section |IT] we 
derive the expansion for pm as a function of m. In Sec- 
tion |III| we illustrate how to estimate entropic exponents 
via fits of with a reweighting of exact enumeration 
data for polymers on square and cubic lattices. In Sec- 
tion |IV| we present Monte Carlo simulations of ELPs con- 
taining a fixed knot and determine the averaged stored 
length in equilibrium as a function of m. The entropic 
exponent 9 is determined for different simple and double 
knots. Finally, in Section |V] we discuss, on the basis of 
the obtained scaling behavior for p,,,,, different possible 
scenarios for knot competitions. 



II. SCALING PROPERTIES OF THE STORED 
LENGTH 

Consider a polymer composed by m monomers with 
lattice coordinates defined by r^, i = 1, 2, ... to. Mul- 
tiple occupancy of neighboring monomers on the same 
lattice site means that we allow configurations for which 
rk = Tk+i = . . . = Tfc+p = s. However if rk-i 7^ s and 
T^k+p+i 7^ s, then no monomers other than those of the 
interval [k,k + p] are allowed to visit the site s. The 
lattice polymer so defined describes a self-avoiding back- 
bone of length < Z < TO. The two extremal cases are all 
monomers occupying the same lattice point (/ = 0) and a 
fully stretched configuration without multiple occupancy 



(Z = to). The equilibrium partition function for an ELP 
with TO monomers is given by 



E 

1=0 



(7) 



u' Zi 



(3) 



where the sum is over the length / of the self-avoiding 
backbone and Zi is the canonical partition function, 
which counts the number of allowed configurations for the 
self-avoiding backbone, and whose asymptotic is given in 
Eq. (Q. The factor (7) in Eq. |3]) counts the number of 
ways the stored length can be distributed over the back- 
bone. For convenience an extra fugacity u per site has 
been added. 

Substituting Eq. ([T]) in (pj) and defining w = jiu, the 
average backbone length (Z) can be computed from 



(0 



d 



(4) 



It is instructive to consider first the case of a partition 
function of the type Zi — /i' in Eq. (jsj), i.e. neglecting 
power-law and correction to scaling terms in Eq. ([T]). In 
this case Eq. ([3| becomes 



Zoo 
m — 



i=0 



J = {1+ujY 



(5) 



Equation (|5| has the following interpretation: the par- 
tition function for a walk of to steps factorizes as each 
monomer can either sit on the backbone (accumulating 
stored length with weight 1) or occupy a free site (with 
average weight w). From Eq. ^ we get the following 
value of the averaged backbone length 



(6) 



We now go back to the full partition function in Eq. (|3| . 
For large m and fixed uj the binomial factor is sharply 
peaked around I = l^o- We approximate the binomial by 
a Gaussian distribution as follows: 



(7) 



J « {i + ujy 



where 



V27rcr2 
muj 

(1+C^)2' 



(7) 



(8) 



The Gaussian approximation differs from the binomial by 
terms which are exponentially small for large to, which 
are of higher order in the large-TO expansion we are in- 
terested in, so they can be safely neglected. We replace 
now the discrete sum in Eq. ([s]) by an integral over all 
lengths, extending the domain of integration in the whole 
real axis: 



7 



{i + ujy 

\/27rcr2 



+ CXJ 



(9) 
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where we have replaced the asymptotic form of Zi as 
given in Eq. ([ij. The replacement of the sum by an 
integral brings corrections in Eq. (|9]) which are of higher 
order in 1/m and for our purposes can be neglected. 

We solve the integral in Eq. ([9| by using a saddle point 
approximation. A simple rescaling I — x loo gives 



z„, - (1 + u^Y' 



(10) 



with 

r(x) 



-{x-\f \0g{x loo) + log (1 + A{X loo)-"^) 



2 mui 
Let X the maximum of r(a;). We have: 



(11) 



27r 



2tt 

7WW\' 



mLj\r"{x)\ 



(12) 



Equation (111 implies that the maximum of T{x) in 



the large-m limit is a; = 1 + C(l/m), giving 

|r"(i)| = 1 



oil 

m 



(13) 



which produces higher-order terms, which we neglect in 
the large-m expansion. In addition: 



mujT{x) = 9 log ^oo + log (l + Al 



OO J 



(14) 



Equations (12 1, (13 1 and (^14| again show that the lead 



ing contribution to the partition function is (l-l-w)™, 
but also that the subleading contribution ^ l^ ^ mP has 
the same entropic exponent 9 of SAWs. From Eq. ([4]) 
we get 



{I) = I. 



1 + 



A 

muj 



A' 



(15) 



and the stored length density ^ becomes 

'1 + 



Pn 



where we defined 



A 



A 

T+A 



Poo = 



1 



(16) 



(17) 



The expansion (16) is valid provided A < 1. The ne- 



glected terms coming from the replacement of the sum 
with an integral, and from the Gaussian integration in 
Eq. (10 1, are of the order 1/m^ (except if A > 1, the 
0{l/rnr) terms would dominate over the l/m^+'^.) The 
value of the exponent 9 can then be obtained from a plot 
of pm vs. 1/m, as the slope in the limit 1/m — > 0. Using 



the high-precision literature values for the connectivity 
constants /z, one obtains a very accurate estimate of Poo- 
With the definition of Poo in ( 17 1 we can rewrite the 
variance (Isl as 



= TOPoo(l - Poo) 



(18) 



This form reveals clearly that the largest a for a given m 
is achieved with poo — 1/2, i.e. with a fugacity u — p^^. 
We can think of this regime as the maximally elastic one. 
In all cases, note that the relative polydispersity a/m of 
the chains goes to zero ~ m~^/^ for m — > oo, hence the 
chain lengths / are narrowly distributed around their av- 
erage (l) . This allows us to use saddle-point approxima- 
tions (see Sect. |v]) and leads to metric properties in the 
universality class of SAWs (e.g. radius of gyration scaling 
as ^ m'^ ^ (O")- Hence, ELPs share both exponents 9 
and v with SAWs. 



III. STORED LENGTH FROM EXACT 
ENUMERATIONS DATA 

As a first illustration of the scaling behavior of 
the stored length as a function of the number of 
monomers m, we consider exact enumeration data for 
SAWs and SAPs on square and cubic lattices, which are 
taken from the published literature |il7j . Enumeration 
techniques provide exact values for the total number of 
SAWs Zi as a function of their length I. We use these 
values for Zi to compute Zm from Eq. (|3|. The stored 
length Pm is obtained from the average (Z)7using Eq. (|2|. 
We have the freedom to choose the value of the fugacity 
u in Eq. (3|. 

Figure 2] shows a plot of pm as a function of 1/m for 
three and two dimensions, obtained by setting u — I 
in Eq. ([3 1. The data converge to the expected asymp- 
totic value, which is poo — 0.175931 (cubic) and poo — 
0.2748643 (square). These are obtained from Eq. ([l7l 
with Lu = up, — p and the following values for the con- 
nectivity constants: p = 4.684044(11) [21? (cubic) and 
p = 2.63815852927(1) (square) L16 



The solid lines in Fig. [2] are the linear terms in the ex- 
pansion of Eq. ( [l6| where the value of 9 is that for open 
walks {9 = 11/32 in d = 2 and 6* ~ 0.157 in d = 3) and 
polygons (6* — —3/2 in d = 2 and 9 ~ —1.76 in d — 3). 
The results show that the linear scaling in 1/m sets in 
already for short polymers (m « 20) . In addition we ob- 
serve that the corrections to the leading scaling behavior 
are stronger for closed walks (empty circles) compared to 
the open walks case (filled circles). 

We performed finite-size extrapolations to obtain es- 
timates of 9 from p,„. The two-dimensional data have 
been extrapolated by means of the Burlisch-Stoer (BST) 
algorithm [32], using the finite-m approximants 



Pm Prn—1 



Pm-i/m - Pm/im - 1) 



(19) 
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FIG. 2. Stored length calculated from exact enumeration data 
for SAWs (bullets) and SAPs (empty circles) on the cubic 
lattice {d = 3) and on the square lattice {d — 2). The solid 
lines are the leading terms in the pm vs. l/m expansion when 
using the expected value for the exponent 6. The asymptotic 
values Poo = 1/(1 + pi) for cubic and square lattices are shown 
as a dashed line. 

which are the ratios between slope and intercept of the 
hne joining the points (l/(m— l),pm-i) and (l/m, p^), 
i.e. they are finite-size estimates of the ratio between paod 
and Poo- The BST algorithm starts with a sequence of N 
elements, and generates iteratively sequences of — 1, 
N — 2 ... elements which are expected to converge faster 
at each iteration step. It involves a free parameter (SI), 
which roughly measures the effective leading correction 
exponent. In our extrapolations, an optimal value of 
was selected requiring a minimal standard deviation of 
the last five sequences generated by the iterative algo- 
rithm. The extrapolations were repeated for different 
values of the fugacity parameter u and the error was es- 
timated from the variation on these values. For three- 
dimensional data, the BST algorithm turned out not to 
be very accurate, particularly for loops. The reason is 
that Pm for small m has some subleading oscillatoric be- 
havior which is not sufficiently damped during the BST 
iterations. The result is that the accuracy of the extrapo- 
lation is poor. For these data we use instead a non-linear 
fit, fixing Poo and keeping 9, A and A as fitting parame- 
ters. 

The extrapolated values for 9 are reported in Table |l] 



TABLE I. Summary of the exponents obtained from the ex- 
trapolation of the approximants drn defined in Eq. ( |19[ ). The 
data are for SAWs and SAPs |17| . The last column gives the 
exact two dimensional data [J. 
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max m [17.1 
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SAW, d = 2 


71 


0.3437(2) 


0.34375(= II) 


SAP, d = 2 


110 


-1.500(1) 


3 
2 


SAW, d = 3 


30 


0.158(2) 




SAP, d = 3 


32 


-1.75(2) 






FIG. 3. (Color online) Examples of knots on the fee lattice: 
(a) a 3i knot (trefoil) with I = 15 steps and (b) a composite 
3i#4i knot with I = 31 steps (the notation fci#fe2 indicates a 
closed polymer ring with two knots, one of type ki and one of 
type fc2). These configurations have been used as backbones 
with stored length m — I for starting the simulations of ELP 
with m > I monomers. 

these are accurate and in good agreement with exact data 
in two dimensions and also with the best numerical esti- 
mates in three dimensions {9 = 0.1573(2) [21] for walks, 
9 = —1.76279(2) for polygons — assuming hyperscaling 
a-2 = -dv with v = 0.587597(7) [Sj), which shows that 
reliable values of the entropic exponents can be extracted 
from the scaling of the stored length. 

IV. ENTROPIC EXPONENTS OF KNOTTED 
POLYMERS 

We now turn to the study of equilibrium properties 
of ELP rings with some fixed topology. Here we will 
show how the knowledge of the stored length p„j can be 
exploited to investigate equilibrium properties of knotted 
polymers. 

We have performed Monte Carlo simulations of ELPs 
on the face-centered-cubic (fee) lattice, with an algo- 
rithm that was recently used to study translocation dy- 
namics [21_ and phase separation in polymer melts [22J. 
The allowed Monte Carlo update moves include repta- 
tion, i.e. the diffusion of stored length along its backbone 
and Rouse-like moves which locally change the backbone 
configuration. (For more details see Ref. |22|). 

The setup of the simulation is as follows. We start 
from a backbone with a minimal number of steps on the 
fee lattice, as those shown in Fig. |3] A total number of 
monomers m are distributed randomly over this back- 
bone. These configurations are then relaxed to equilib- 
rium. Typically m is much larger than the initial length 
(we simulated polymers with m up to 2000) so that re- 
laxation to equilibrium corresponds to an expansion of 
the backbone. The Monte Carlo moves preserve the knot 
topology imposed initially. Once equilibrium is reached 
we start the sampling of the stored length density pm- 

An additional weight (equal to 4) is introduced for 
moves that accumulate monomers on the same lattice 
point, which corresponds to a fugacity factor u = 1/4 in 
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FIG. 4. (Color online) Plot of the average equilibrium stored 
length density pm as a function of the inverse monomer 
number 1/m for closed polymers with some fixed topology. 
The simulations were extended to polymers of lengths up to 
m = 2000. From top to bottom the data refers to: unknotted 
ring, 3i, 4i, 5i and 62 knots. The two bottom data set cor- 
respond to configuration with two knots: 3i#3i and 3i#4i, 
respectively. The horizontal dashed line is poo ~ 0.28498, as 
expected from Eq. (201. Inset: zoom of the asymptotic re- 
gion. Straight lines represent pooi^ — 0/m): the dotted line 
corresponds to the conjectured value 9 — —Sv, the dashed 
line to ^ = —Sf + 1, and the dot-dashed line to 9 = + 2. 



Eq. ([3]) . This leads to the following asymptotic value for 
the stored length density: 



1 + /i/4 



0.28498(1) 



(20) 



(20 



This value is shown 
b ig. |4] We note that 



where the numerical value is obtained by considering the 
most accurate available estimate n = 10.0362(6) [9 for 
the connectivity constant of SAWs on fee lattice [? ]. 

Figure |4] shows the scaling behavior of pm as a func- 
tion of 1 /m for an unknotted polymer ring, for single and 
double knots. All data converge asymptotically to the 
value Poo obtained from Eq 
as a dashed horizontal line in 

the approach to poo of the numerical data for unknot 
ted rings is quite different for those of rings with knots 
(a detail of the asymptotic region is shown in the inset 
of Fig. [4]): the data for the unknotted topology approach 
the asymptotic value with a clear 1/m scaling behavior. 
For topologies with knots instead there is a pronounced 
curvature in the pm vs. 1/m plot, deriving from strong 
corrections to scaling. These corrections are stronger for 
an increasing knot complexity and for an increasing num- 
ber of knots. The shortest length ^i„in of a knot on a 
lattice is a good indicator of its complexity, and in this 
model for m = I — ^min by definition the chain can only 
be fully stretched, i.e. Pmrnm — 0. Our data show that 
the crossover from this initial topological stretching to 
the asymptotic regime ^ m^^ grows quickly with the 



smaller than that of the simplest knot (the trefoil with 
"T-min = 15) explains why corrections to scaling are neg- 
ligible for unknotted chains. 

We estimated the entropic exponent 6 using the scal- 
ing behavior predicted by Eq. (16 1. In the unknotted 



case due to the manifest absence of curvature of the 
data, we restricted ourselves to a linear fit setting A = 
and Poo = 0.28498 in Eq. (16) and using 9 as the only 
free parameter. The fit, restricted to m > 200, yields 
9 = —1.76(3), confirming that the entropic exponent for 
rings with fixed unknotted topology is identical to that 
for SAPs with no topological constraints. 

A closer look at the data reveals that the stored length 
density for knots 3i, 4i and 5i is non-monotonic. As the 
data asymptotically approach poo from above, Eq. ( 16 1 



implies a negative value of the exponent 9. We performed 
a non- linear three-parameters fit to the data based on 
Eq. (16 1: 9, A and A are fitting parameters while we fix 
Poo = 0.28498, as predicted by Eq. (|20|. The results of 
the non-linear fits are given in Table |II| The estimated 
exponent changes sign from single knot (9 < 0) to double 
knots (9 > 0). A range of correction-to-scaling exponents 
A providing optimal fits were selected and these are given 
in the third column of Table |TT1 Error estimates for 9 
refiect the variability in 9 from the different values of 
A used in the analysis. For the knots studied, the most 
accurate estimate for 9 is that of the 3i knot, yielding 9 = 
—0.75(5). The error increases with the knot complexity. 
For single knots we also note a change in the range of 
correction-to-scaling exponents from A « 0.6 for the 3i 
knot to A w 1.1 in the other knots. 

It has been suggested |13 | IT4 j l35| that for a knot k with 
TTfc prime components, the entropic exponent is given by 
^fc = ^ + ""fe ! where 9 is the exponent for a polymer ring 
without fixed topology. If this is the case we expect for 
a single knot an exponent 9 = —3^+ 1 = —0.76 while for 
double knots 9 — —iv + 2 = 0.24 (these conjectured val- 
ues are shown as dashed lines in the inset of Fig. |4|. The 
idea behind this suggested scaling is that localized knots 
are like sliding entities, which can occupy any of the I 
sites of a chain, thus contributing entropically with a fac- 



TABLE II. Summary of the estimated entropic exponents ob- 
tained from the scaling behavior of the stored length with a 
three parameters fit (5, A and A in Eq. ( 16 1). The asymptotic 
value Poo is kept fixed. The last column shows the range of 
polymer sizes used in the fit. 



value of TTln 



In this view, the fact that for unknot- 
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Knot type 
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A 


range of m 


unknotted 


-1.76(3) 




> 200 


3i 


-0.75(5) 


0.5-0.7 


> 400 


4i 


-0.6(1) 


0.9-1.2 


> 300 


5i 


-0.5(2) 


0.9-1.2 


> 300 


52 


-0.3(3) 


0.9-1.2 


> 300 


3i#3i 


0.4(2) 


1.2-1.4 


> 300 


3i#4i 


0.8(2) 


1.2-1.4 


> 300 
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FIG. 5. (Color online) Plot of fm.m^'^^ vs. logjg m, with 
A = 1/2. The data tend to a constant for large m, which, as 
discussed in the text, is consistent with a correction-to-scaling 
exponent of A « 0.5. 



tor I in the partition function. Our numerical results fully 
support this conjecture for the single 3i knot and also for 
the double 3i#3i knot. Results for the other knots seem 
to overestimate with respect to the conjectured values. 
It is likely that the deviations from the conjectured val- 
ues are due to strong finite-size effects. An indication of 
this is the value of the correction-to-scaling exponent ob- 
tained from the fits, which, with the exception of the 3i 
knot, is estimated as A « 1. Renormalization group ar- 
guments IJ4 for magnetic 0{N) models, which map into 
polymer models in the limit [2|, predict instead 

A « 0.55, and Clisby [5 finds A w 0.528(12) in simula- 
tions of very long SAWs (this is in agreement with the 
range of values obtained in the extrapolations of the nu- 
merical data for the 3i knot, see Table [n]). We also re- 
mark that a value A > 1 is at odds with the expansion 
of the stored length of Eq. ( 16 1 in which it was implicitly 
assumed A < 1, the term would be otherwise 

dominated by 0(l/m^) corrections, which were neglected 
in the computation of pm leading to Eq. (16 1. 



In the sequel we fix the entropic exponents to the con- 
jectured values = + TTk and subtract from pm the 
constant and leading correction in 1/m as 



fn 



1 



m 



(21) 



swollen polymers [3|. In addition the constant A is nega- 
tive and its magnitude quickly increases with knot com- 
plexity. This is also visible in Fig. |4] as the effect of in- 
creasing A is that of producing an increased curvature in 
a plot of pm vs. 1/m. It is perhaps not surprising that 
finite-size effects increase with the knot complexity, as 
more complex knots are expected to occupy a larger por- 
tion of the polymer. In table |ITT] we list our estimates of 
A and B, obtained by means of linear fits to data in the 
form fm^n^'^ vs. m"*^'^. The values of B are almost two 
orders of magnitude larger than those of A, explaining 
the fitted (effective) leading exponent A « 1. 



V. KNOTS COMPETITION 

In this Section we discuss entropic competition be- 
tween knotted polymers in the context of ELPs. The 
idea of entropic competition between polymers with var- 
ious constraints was introduced in Ref. [25j as a direct 
way to estimate polymer entropic exponents from canon- 
ical simulations. This idea is sketched in Fig. |6] and 
can be implemented in various ways. One can consider, 
for instance, a polymer loop divided in two sides by a 
wall (Fig. [6f^a)); the two sides exchange monomers via 
sufficiently small holes such that the knots cannot pass 
through. The exchange can also occur through a ficti- 
tious "wormhole" ^25j, as shown in Fig. [6]jb). The poly- 
mers at the two sides of the wall or those exchanging 
monomers through the wormhole do not interact with 
each other. When exchanging monomers the length of 
each loop fiuctuates, while the total length is fixed to a 
constant L. The method [25_ is based on the analysis of 
the equilibrium distribution of lengths of the two sides. 
For ordinary polymers one expects that the length / of 
one polymer ring is distributed according to 



p{l) ^ Z\'>Z^ 



(1) 7(2) 



(23) 



where the two Z's are the loop partition functions given 
in Eq. ( |l|). The main point is that the dependence on p 
in Eq. ( |23| is irrelevant as L is fixed, whereas from the 
analysis of the shape of the probability p{l) as a function 
of I it is possible to fit the values of the entropic exponents 
6*1 and 62 of the two loops ^25j. 



For this quantity we expect the following scaling behavior TABLE III. Fits of A and B, assuming A = 1/2 in Eq. (22i 



m 



A 

l+A 



B 



(22) 



where a next-order l/m^ term has been added. 

Figure [5] plots fm'm}'^^, where we set A = 0.5, as a 
function of m. The fact that this quantity approaches a 
constant value for large m supports an estimate of the 
correction-to-scaling exponent A w 0.5, as expected for 



Knot type 


A 


B 


3i 


-0.78 


-33 


4i 


-1.9 


-63 


5i 


-3.5 


-73 


52 


-4.1 


-80 


3i#3i 


-1.4 


-130 


3i#4i 


-3.3 


-173 



(b) 



FIG. 6. (Color online) (a) Sketch of an entropic competition: a portion of the ring polymer, with two knots, is constrained to 
stay on the left half-space (holes are small enough to forbid more than one monomer at a time to pass), the remaining part 
has one knot and is on the other side of the wall. The total length A'^ of the chain is constant but the lengths m and N — m oi 
the two subchains can fluctuate, (b) The virtual version (without the wall) of the same competition: the two polymers swim 
in separate dilute solutions and are coupled via a "wormhole" trough which they can exchange a monomer (hence not a knot) 
at a time. 



Entropic competition without a wall 



We first consider the case depicted in Fig.|6|b), and we 
discuss a few representative examples. If the two loops 
both have negative entropic exponents (^i , ^2 < 0) , then 
one expects a,p{l) as depicted in Fig.|7|^a) (thick line and 
shaded area), whereas the case 01,6*2 > is depicted in 
Fig.jrjjb) (same notation; in these figures, for convenience 
we show the distribution p{l/L), which is just a rescaling 
of p{l)). The thin lines in Fig.lWa) and (b) show sketches 
of finite-L distributions of p(/) for increasing L: partic- 
ularly interesting is the scenario depicted in Fig. |7|^a), 
which shows a drastic change of the shape of the dis- 
tribution from a finite L to the the limit L —> 00. We 
will discuss here how some of these features can be un- 
derstood from the analysis of the stored length densities 
/5m'' and pm'' of the competing loops. 

Let us consider p(m), the probability of finding m 
monomers in one of the two entropically competing ELF 
loops. This quantity scales as 



p{m) 



(24) 



where the two Zs are the partition functions of the two 
competing ELPs at fixed monomer numbers m and with 
fluctuating lengths. To find the most probable value of 
the monomer number m* observed in the entropic com- 
petition setup, we maximize the entropy 



Sra--kB l0gZ« 



kB log4'-r 



(25) 



(kB is the Boltzmann constant). The partition functions 
of ELPs in Eq. ([3| are expressed as a sum over all lengths 
< m < Z. The sum is however dominated by a charac- 
teristic value of l*{m) obtained from the condition 



dZ„ 



dl 







(26) 



Dh 0- 




(b) 








short L 










L — >oo 





U.4 0.6 

l/L 



FIG. 7. (Color online) Examples of probability distributions 
of loop lengths for two polymer loops exchanging monomers 
as in Fig. |6f^b) , in the case of negative (a) and positive (b) 
entropic exponent 9 for both loops. In (a) the competition is 
between two 3i knots, the thick line (boundary of the shaded 
area) is the distribution for the limit of long L while the other 
ones are for two short L's. In (b) the competition is between 
a 3i#3i knot and a 3i#4i knot, with the same notation. 



where we defined 



(7) 



u'Z, 



(27) 



Now assuming that Zm is dominated by a single value of 
r{m) we can compute the total derivative in ni of log Zm 



d log Z„ 
dm 



'^log^m,C(mO _ (^l0g^m,;*(m,) 

m! 



dm 



dm 



dm ^ (m — l*{m))\ 



logpm. (28) 



In this derivation we used Eq. ( 26 1 , so in the total deriva- 



tive with respect to m we can ignore the m-dependence 



1=1' {m) 



coming from l*{m). Combining Eqs. (25 1 and (28 1 we 
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FIG. 8. (Color online) Plot of stored-length densities vs. m 
and N — m in the entropic competition setup, for four different 
knotted chains: (a) 3i#3i vs. 3i#4i, (b) unknot vs. unknot, 
(c) 3i vs. 3i, and (d) unknot vs. 3i#3i. The intersection 
points of the densities are highlighted by arrows and cor- 
respond to local maxima (filled arrows) or minima (empty 
arrows) of the total entropy of the two competing loops. Hor- 
izontal dotted lines indicate poo- 



find that the extremum of the entropy Sm of the com- 
peting rings is given by the value of m* for which 



(1) 



(2) 



(29) 



To find TO* one can plot p'm and p)^'_„i vs. to and N ~m 
in the same graph: each intersection point between the 
two curves is an extremum of Sm ■ To decide whether this 
is a local maximum or minimum one analyzes the second 
derivative 



,(2) 



dw? 



1 



(1) 



dp 
dm 



^ (2) 

d-PN-m 

d{N - m) 



(30) 



In Fig. [8] we show some plots of the stored-length den- 
sities for the two competing loops containing knots. The 
two loops have a total number of monomers equal to 
N = 2000 and the data are those shown in Fig. [4] but 
now plotted as function of to and N — m. As seen in the 
previous section, the stored length density can be non- 
monotonic in TO for some knotted configurations, which 
can produce various scenarios where up to three intersec- 
tion points are possible. 

Figure [Sjja) shows the example of two competing dou- 
ble knots. In this case there is a single intersection point 
and the analysis of the first derivatives of p shows that 
this point is a local maximum for the entropy (Eq. ( |30[ )). 
The probability distribution of monomers (or lengths in 
the canonical setup) will have a single maximum at some 
intermediate to* , as shown in the example of Fig. [7jb) . 
In the case of two unknotted loops (Fig. [sj^b)) the in- 
tersection point TO* is a minimum for the entropy, hence 
the probability distribution for m will be maximal at the 



edges and minimal at to*, as for the thick line in Fig.|7|^a). 
The most interesting case is that of competition between 
loops with non-monotonic p's. This case is illustrated 
in Fig. [sjc). The three intersection points are a cen- 
tral local minimum of the entropy enclosed by two local 
maxima. The probability distribution of lengths is like 
that depicted as a dashed line in Fig. |7|^a). It is easy to 
see that if the total number of monomers decreases (this 
corresponds to shift one of the two p's along the horizon- 
tal axis) there will be only one intersection point. This 
generates a probability distribution with a single maxi- 
mum for TO* (thin dense line in Fig.|7ja)). Interestingly, 
the length distributions obtained from Monte Carlo si- 
mulations (25] of competing off-lattice flexible rings with 
simple knots give, for sizes up to 200 monomers, concave 
distributions, contrary to the expectations of negative 0's 
from the conjecture of Ref. [IS] , which would instead cor- 
respond to a convex (i.e. with a minimum in the middle) 
shape. The non-monotonicity in to of the stored-length 
density pm explains this drastic change in behavior in 
finite-size data. 

To complete the discussion we consider next an exam- 
ple where no intersection point is present (Fig. [sjd)). In 
this case one has to resort to the full form of Sm'- from 
the scaling of partition functions of SAPs, the probability 
of a state with m monomers on the side with no knots 
is expected to scale as m~^'^ ~ to"-'^-^^ with a cutoff at 
m < N . It implies that the average length of the unknot- 
ted subchain (to)^ ~ N'^ ~ JV"-^'^ is weakly scaling with 
N , and at least in this case the competition is clearly 
in favor of the side with knots. This reminds us that 
the full statistics given by Sm would often be necessary 
to compute average quantities, and that the maxima are 
only indicative elements. Nevertheless, we have seen that 
the density of stored length is a useful quantity for under- 
standing the basic properties of the entropy of competing 
knotted chains. In particular, knowing it and its short- 
features helps to interpret the numerical results and to 
distinguish preasymptotic scalings from asymptotic ones. 



2. Entropic competition with a wall 

Let us finally go back and reconsider briefly the en- 
tropic competition of knots divided by a wall, as in 
Fig. [6]^a). The main difference is that the basic expo- 



nent of the unknotted chain should be 



— —dv 



The additional index is connected to the constraint 
of having a monomer of a loop confined close to a hard 
surface. The formula is an application of Duplantier's 
general theory of polymer networks |43[ HI) . We use this 
theory also to extract a'2 from the data in Ref. (45], ob- 
taining ~ —0.95. This means that Og « —2.71. Again 
the full zoology of possible competitions could be simply 
discussed by repeating the above reasoning, once data of 
p„i for ELPs close to a wall are generated. We reserve 
this investigation for a future work. Let us just note that 
the condition Og+Tik > 0, associated with a single maxi- 
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mum of the entropy Sm at 1 <C m* ^ iV, is now met for 
a minimal number of prime knots tt^. = 3 per loop, i.e. 
one more than we needed in the case without the wall. 
Thus, the wall separating the chains has somewhat the 
effect of repelling entropically also the knots. 



VI. CONCLUSIONS 

In this paper we studied the scaling properties of a 
class of polymers, which we have referred to as elastic 
lattice polymers (ELPs). These polymers can accumu- 
late stored length along their backbone, by lifting the 
self-avoidance condition for neighboring monomers. The 
length I of their backbone fluctuates, whereas the total 
number of monomers m is fixed. Differently from true 
grand canonical polymers, however, the backbone length 
is bounded to Z < m, and fluctuates around (/) with fluc- 
tuations ^ ^/m. 

ELPs were used in the past to study the dynamics of 
polymer melts [22 and pore translocation dynamics |21| . 
but their equilibrium behavior has received little atten- 
tion. In this work, we used ELPs to investigate entropic 
exponents of knotted polymer rings. The calculation of 
polymer entropic exponents from classical Monte Carlo 
simulations of canonical self-avoiding rings is rather cum- 
bersome: one has either to employ complex grand canon- 
ical sampling, or to resort to the so-called atmospheres 
method [19 or entropic competition methods [2S" . ELPs 
are well suited to this type of problem. Firstly, the under- 
lying Monte Carlo dynamics can be based solely on local 
moves (thus conserving the knot topology) and the pos- 
sibility of accumulating length along the backbone facil- 
itates the sampling of different conflgurations compared 
to canonical self-avoiding rings. The ELP explores new 
conflgurations through sliding moves along the backbone. 
Secondly, we have shown that there is a natural variable 
associated to ELPs which is the stored length density 
Pm (see Eq. [2|, which measures the average fraction of 
monomers accumulated on the backbone. We have de- 
rived an expansion for pm in the limit m — >■ cxd, where 
Pm converges to a value that depends on the connectiv- 
ity constant of the ordinary lattice polymers. The next 
leading behavior is of the order 6/m, with 9 the entropic 
exponent of the polymers. This allows one to estimate 



entropic exponents from the scaling analysis of pm- As 
examples of application of this, we estimated entropic ex- 
ponents of swollen polymers in d = 2 and = 3, and of 
polymers with various types of knots. Comparing the re- 
sults with conjectured values of these exponents, we flnd 
a clear agreement at least for the simplest knot studied. 
For more complex knots the agreement is only marginal, 
due to flnite-size effects quickly increasing with the knot 
complexity. 

One of the advantages of the stored length analysis 
is that correction-to-scaling effects are directly visible in 
Pm vs. 1/m plots as they appear as deviations from a 
linear scaling behavior. Our analysis showed that flnite- 
size effects become stronger with the knot complexity 
and with the number of knots. Similar result have been 
observed by Janse van Rensburg and Rechnitzer [19J. 
These authors estimated the connectivity constant and 
entropic exponent of lattice polymers via the atmosphere 
method [18j , where, roughly speaking, atmospheres are 
the loci where the polymer can be expanded and con- 
tracted. Interestingly, there is a similarity between the 
scaling of the average atmospheres and that of the stored- 
length density of ELPs discussed in this paper. 

With simulations of ELPs we have shown how impor- 
tant are corrections to scaling in the statistics of knotted 
polymers: their equilibrium properties in entropic com- 
petition can be understood from coexistence diagrams of 
stored lengths of ELPs. The non-monotonicity of the 
stored length density as function of l/L explains some 
features of the competing rings observed in canonical 
Monte Carlo simulations [2S] which were poorly under- 
stood before. 

Summarizing, the elastic lattice polymer is a simple 
model sharing critical exponents with the self-avoiding 
walk, but it has an additional "elastic" degree of freedom 
in its fluctuating length, which offers numerical advan- 
tages and additional theoretical tools to derive critical 
exponents of polymers. Thus, the ELP is a valid alter- 
native to classical lattice models for studies in polymer 
physics. 
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